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Abstract 

LEED experiments show that Li adsorbed at Cu(lOO) surfaces at 
room temperature induces a (2x1) missing row substrate reconstruc- 
tion while adsorption at lower temperatures, T=180 K, results in an 
unreconstructed Cu(100)+c(2x2)-Li overlayer structure. Substrate 
reconstruction has not been observed for Na nor for K adsorption. 
In order to study the specific reconstruction behavior of the Li adsor- 
bate ab initio DFT calculations have been performed on Cu(100)+Ad, 
Ad = Li, Na, K systems at coverages 0Ad = 0.25 — 0.5 with and with- 
out reconstruction. The calculations show that the (2x1) MR recon- 
structed surface lies energetically above the ideal (1x1) surface by 
0.2 eV per unit cell. However, alkali binding is stronger in the MR 
geometry as compared to that of the ideal surface where the increase 
in bond strength becomes smaller in going from Li to Na to K. As 
a result, the MR reconstructed and the overlayer adsorbate systems 
are energetically very close for Cu(100)+Li while for Na and K the 
overlayer geometry is always favored. 

1. Introduction 



It is well known from experiment that alkali adsorption at fcc(llO) surfaces 
induces a (2x1) missing row (MR) substrate reconstruction (for reviews see 
[]T], 0). Recently, this type of adsorbate induced reconstruction was also found 
on fee and bcc (100) surfaces. Examples are Au(100)+K |^, Au(100)+Na 
i, Ag(100)+K |], and Cu(100)+Li [|, 0. Here the Cu(100)+Li system is of 
particular interest since it exhibits different reconstruction behavior depend- 
ing on the adsorption temperature H, 0. At low temperatures, T < ISOi^', 
a c(2x2) overlayer (OL) structure is observed for 6Li ~ 0.5 without indi- 
cations of a substrate reconstruction. However, subsequent heating of the 
adsorbate system to room temperature leads to a (2xl)MR reconstruction 
of the Cu(lOO) substrate which is conserved even if the system is cooled back 
to 180 K. This type of adsorbate induced substrate reconstruction seems 
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unique for Li adsorption on Cu(lOO) and has not been found for larger al- 
kali adsorbates such as Cu(100)+Na g], Cu(100)+K i, 0, or Cu(100)+Cs 



T0| . A theoretical understanding of this behavior is still lacking. In the 
present theoretical study we have studied the underlying mechanisms by ex- 
amining differences in the alkali adsorption on ideal (1x1) and (2xl)MR 
reconstructed Cu(lOO) in order to explain the specific behavior found for Li 
adsorption. 



2. Computational Details 

Geometric structures of the Cu(lOO) substrate and of the Cu(100)-|-Ad, 
Ad = Li, Na, K adsorbate systems are described in the repeated slab ge- 
ometry. The Cu substrate, in its unreconstructed (1x1) and (2xl)MR re- 
constructed form (see Fig. 1), is approximated by 7 and 9 layer slabs which 
are separated by vacuum corresponding to 5 substrate layers. This has been 



found to yield a reliable surface representation ||Tl|. The alkali adsorbate 
layers (Li, Na, K) are added at the top and bottom of each substrate slab 
where for coverages = 0.25 a lateral (2x2) and for = 0.5 a c(2x2) and 
(2x1) supercell geometry is assumed as illustrated in Fig. 2a-f. The geome- 
tries of the surface systems are further optimized by allowing all atoms in the 
unit cells to relax according to the force field obtained from the calculations 
described below. 

The electronic structure of the surface systems is calculated within the 
density functional theory (DFT) formalism using the local density approxi- 
mation (LDA) [l^, [l^ as well as gradient corrected functionals (GGA-II) [l^ 
for exchange and correlation. Total energies and derived quantities are ob- 
tained with the full-potential linearized augmented plane wave (FP-LAPW) 
method [|T5| implemented in the WIEN93 program |1^ where 45 k-points in- 



side the irreducible part of the Brillouin zone are used and the energy cut-off 
for plane wave expansions in the interstitial region is set to 9.58 Ry. Ge- 
ometry optimizations of the surface systems are based on calculations of 
forces acting on corresponding nuclei [|l7l where equilibrium was assumed 



for forces below 3 mRy/bohr. A comparison of the geometric and energetic 
results from calculations using the LDA and the GGA-II approach shows 
only small differences which do not affect the present conclusions. Thus, we 
restrict ourselves in the following to results from our LDA studies. 
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3. Results and Discussion 

The (2xl)MR reconstructed Cu(lOO) surface, shown in Fig. 1, is created 
from the ideal unreconstructed surface by removing every second row of sur- 
face atoms. As a result the top-most atoms of the reconstructed surface 
experience a different electronic environment which affects the surface ge- 
ometry. This becomes obvious from Table la which gives numerical results 
of the present geometry optimizations. For the unreconstructed Cu(lOO) 
surface the interlayer spacing between the first and second surface layer, 
di2 (= 1.76 A in the ideal bulk geometry), is decreased, A(ii2 = —3.2%, and 
that between the second and third layer, 0^23, slightly increased, A(i23 = 0.2%, 
in good agreement with experimental |]18[ and other theoretical [|19| results. 



In contrast, the (2xl)MR reconstruction of pure Cu(lOO) yields much larger 
changes in the interlayer spacings: di2 is decreased by almost 10% and ^23 
by 2%. In the present calculations the surface energy 7 is defined by 

7 = [Etot(slab, n) - ^tot(bulk, n)]/2 (1) 

where the total energy (per unit cell), ii^tot(slab, n), includes surface relax- 
ation and reconstruction of the slab of n layers while £'tot(bulk, n) is obtained 
for the corresponding slab of the ideal bulk. For ideal (1x1) Cu(lOO) the cal- 
culations yield 7 = 0.8 eV while for the (2xl)MR reconstructed surface one 
finds 7 = 1.0 eV per (1x1) unit cell. Thus, the MR reconstruction requires 
0.2 eV per unit cell for Cu(lOO) which is about 10 times larger than the 



energy needed for the MR reconstruction of the Cu(llO) surface pO 



Table la compares results from geometry optimizations on the different 
Cu(100)+Ad, Ad = Li, Na, K adsorbate systems. Here the left part refers to 
overlayer structures on unreconstructed Cu(lOO), (2x2) for adsorbate cover- 
age e = 0.25 (Fig. 2b) and c(2x2) for 9 = 0.5 (Fig. 2c) The right part shows 
results for adsorption on (2xl)MR reconstructed Cu(lOO), with (2x2) adsor- 
bate geometry for 6 = 0.25 (Fig. 2e) and with (2x1) for 6 = 0.5 (Fig. 2f). 
Test calculations on Cu(100)+Li with different lateral adsorbate sites showed 
that sites of maximum coordination close to Cu substitutional sites were al- 
ways energetically preferred in accordance with Figs. 2b, c, e, f. Therefore, 
this lateral geometry was also assumed for the Cu(100)+Na and Cu(100)+K 
systems. 

For Cu(100)+Li the overlayer geometry (Figs. 2b, c) yields stable Li cen- 
ters at -Zcu-Li = 1.8 A above the substrate surface with only small differences 
between coverages Gli = 0.25, 0.5. This is in good agreement with recent 
experimental LEED results, cp. Table lb. Further, substrate relaxation. 
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quantified by A(ii2, A(i23, seems to be little influenced by adsorption as ev- 
idenced from a comparison of the theoretical Cu(100)+Li with the Cu(lOO) 
data. The overlayer geometry results for Cu(100)+Na and Cu(100)+K are 
found to be rather similar to those of Cu(100)+Li. Here the adsorbates sta- 
bilize for both coverages at distances 2;cu-Na = 2.1 A and ^cu-k = 2.5 A 
above the unreconstructed Cu(lOO) surface. The distance values are larger 
than 2;cu-Li given above which is consistent with the increased atom size in 
going from Li to Na to K. 



Geometry differences between the different alkalis are much more pro- 
nounced in the results for adsorption on the (2xl)MR reconstructed Cu(lOO) 
surface. For Cu(100)-|-Li the Li atoms arrange along the Cu troughs (see 
Fig. 2e, f) stabilizing only slightly above the first Cu substrate layer, 
zcn-u = 0.4 A for Qu = 0.25, ^cu-Li = 0.2 - 0.3 A for Qu = 0.5. The two 
different zqu-u values given in Table la for Qu = 0.5 refer to results of ad- 
jacent Li centers in rows along the Cu troughs indicating a buckling of the 
adsorbate layer. This is reasonable due to the fact that the Cu-Cu nearest 
neighbor distance which determines the lateral distance of the adsorbate rows 
is smaller than typical Li metal nearest neighbor distances. The increased 
substrate relaxation calculated for the pure (2xl)MR reconstructed Cu(lOO) 
surface, see above, exists also in the Cu(100)+Li adsorbate system with vari- 
ations of 3%, see Ad^, Ad23 of Table la. A comparison of the theoretical 
relaxation data with those of the LEED experiment on Cu(100)+Li |TD 



cp. Table lb, yields good agreement. The geometries calculated for Na and 
K adsorption on (2xl)MR reconstructed Cu(lOO) differ considerably from 
those of the Cu(100)-|-Li system. For the lower coverage, 9 = 0.25, the ad- 
sorbates stabilize at distances 2;cu-Na = 0.86 A and zqu-k = 1-7 A above the 
substrate surface which is larger than the -Zcu-Li value and explained by the 
different atom sizes. For = 0.5 the adsorbate layer buckling is much more 
pronounced for the larger alkalis. While in Cu(100)+Li the hight difference 
between adjacent adsorbate centers is only Azcn-u = 0.11 A the calculations 
yield A^cu-Na = 0.31 A and A^cu-k = 2.31 A. In fact, in Cu(100)+K the 
theoretical buckling is so large that the concept of a smooth adsorbate layer 
becomes questionable. 



Fig. 3 shows color coded contour plots of the charge rearrangement in the 
Cu(100)+alkali systems due to adsorption on the (2xl)MR reconstructed 
substrate. The rearrangement characterizes the electronic coupling of the 
adsorbates with the substrate. It is quantified by the valence electron differ- 
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ence between the adsorbate system and its separate components 

Apw(r) = Pva/(Cu(100) + Ad, r) - pw(Cu(100), r) - pw(Ad, r) (2) 

where positive Apva/ values indicate charge accumulation and negative val- 
ues depletion as a result of the adsorbate-substrate interaction. The plots 
refer to Li (Fig. 3a), Na (Fig. 3b), and K (Fig. 3c) coverages 9 = 0.5 and 
show Apva«(i^) for ^ planar cut perpendicular to the surface along the sub- 
strate troughs with Cu centers marked by circles and adsorbate centers by 
squares. Obviously, the charge rearrangement leads, for all three systems, 
to a depletion of charge above the (buckled) adsorbate layer which is con- 
nected with the overall positive charging of the adsorbates at the surface. 
The latter has been confirmed by charge integration and is to be expected 
from simple electronegativity arguments. Further, the plots show an almost 
homogeneous charge accumulation in the region between adsorbate and sub- 
strate suggesting no directional adsorbate-substrate bond formation. Thus, 
the alkali+Cu interaction is determined to a major degree by electrostatic 
and metallic contributions. Corresponding Apvai(r) plots for the adsorption 
on unreconstructed Cu(lOO) lead to identical conclusions. These results con- 
firm the complexity of alkali adsorption on metal surfaces which has been 



discussed in the literature for a long time. For reviews see Refs. [pT] , p2 . 



The relative stability of the Cu(100)-alkali systems for adsorption on the 
(2xl)MR reconstructed vs. the unreconstructed Cu(lOO) substrate can be 
determined from total energy studies. For this purpose we define a recon- 
struction energy Es{Q) by 

Es{Q) = 7rec - 7unrec + -Drec(0) " -Dunroc(0) (3) 

where 7rec, 7unrec dcuotc surface energies of the pure substrate with and with- 
out reconstruction while Drec(0), -Dunrec(0) are adsorbate binding energies for 
the two substrate geometries obtained from respective total energies. Thus, 
Es{Q) includes the energy required to reconstruct the substrate as well as 
the difference of adsorbate-binding with and without reconstruction. Positive 
Es{Q) values suggest that the unreconstructed adsorbate system is energet- 
ically preferred while negative values yield the reconstructed system to be 
favored. Fig. 4 shows the result of the present calculations in a level diagram 
of Es{Q) for adsorbate coverages = 0,0.25,0.5 (0 = denotes the pure 
substrate). On pure Cu(lOO) the reconstruction energy amounts to Eg = 0.2 
eV per unit cell. This energy is reduced to 0.03 eV for Li and Na and to 
0.05 eV for K adsorption at coverages = 0.25 which demonstrates that the 
adsorbates are always bound more strongly to the (2xl)MR reconstructed 
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than to the unreconstructed Cu(lOO) substrate. At an increased coverage 
G = 0.5 the reconstruction energy is further decreased to 0.02 eV for Li ad- 
sorption. In fact, a conservative estimate of numerical errors in the calcula- 
tions, yielding ±0.02 eV, suggests that for Cu(100)+Li the reconstructed and 
unreconstructed systems become equally stable {Eg ~ eV) for 9 = 0.5. In 
contrast, the reconstruction energy leads to positive values for the two larger 
alkah adsorbates, ^,(6 = 0.5) = 0.14 eV for Na and ^,(9 = 0.5) = 0.19 eV 
for K. This indicates that for Cu(100)-|-Na and Cu(100)+K adsorption on 
the unreconstructed substrate is energetically preferred which confirms the 
experimental findings ^ 

4. Conclusions 

The present DFT total energy studies reveal interesting details of the adsor- 
bate induced substrate reconstruction in the Cu (1 00 )-|- alkali systems. Ge- 
ometry optimizations show that in Cu(100)+Li the adsorbate can stabi- 
lize on the (2x1) MR reconstructed surface above the substrate troughs 
where the adsorbate-substrate binding is stronger compared to adsorption 
on the unreconstructed surface. As a result, the adsorbate systems with 
and without substrate reconstruction become equally stable for coverages 
G = 0.25 — 0.5 and transitions between the two becomes possible with vary- 
ing temperature confirming the experimental findings |^. In contrast, 
calculations on the Cu(100)-|-Na and Cu(100)+K systems show that, for 
coverages 6 = 0.25 — 0.5, adsorption on the unreconstructed Cu(lOO) sur- 
face is always energetically preferred over that on a (fictitious) (2x1) MR 
reconstructed surface, again confirming the experiment [|, ^ |10|- Based on 
the calculations, this different behavior between the alkali adsorbates can 
be explained in parts by simple geometric effects. The larger alkalis, Na, K, 
cannot stabilize close enough to the reconstructed substrate surface and their 
atom sizes lead to strong layer buckling for higher adsorbate coverages which 
makes adsorption on the reconstructed substrate very unlikely to occur. It 
should be emphasized that the present calculations refer to static equilibrium 
structures of Cu(100)+alkali which can be compared for different constraints, 
overlayer or missing row. Transition states between the different structures 
and respective energy barriers require the knowledge of appropriate reaction 
paths which have to be assumed or obtained from separate ab initio molecular 
dynamics calculations. 
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Tables 



Table 1. 

Geometry results of the Cu(lOO) substrate and of the Cu(100)+Ad, 
Ad = Li, Na, K adsorbate systems with and without substrate reconstruc- 
tion and for adsorbate coverages B = 0.25 — 0.5. Table la lists theoretical 
data from the present LDA calculations and from those of ref. 
Table lb compares experimental data of refs. [0 
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while 

^ ^ ^ ^ ^ ^. For 

definitions of the (differential) perpendicular distances, Zc^_p^d, Adu, ^^(^23, 
see text. 
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Figure Captions 

Fig. 1 

Perspective view of the Cu(lOO) surface showing the ideal (1x1) and the 
(2x1) missing row reconstructed geometry. The two-dimensional elemen- 
tary cells are marked by rectangles for each geometry. 

Fig. 2 

Top view of the pure Cu(lOO) and of the adsorbate covered Cu(100)+alkali 
surfaces illustrating different surface geometries considered in the present 
study. The Cu substrate atoms are shown as dark (1st layer) and hght (2nd, 
3rd layer) gray balls while the alkali adsorbate atoms are indicated by red 
balls. The upper part of Fig. 2 refers to adsorption on the unreconstructed 
Cu(lOO) surface showing (a) the ideal (1x1) substrate, (b) (2x2) adsorption 
(© = 0.25), and (c) c(2x2) adsorption (0 = 0.5). The lower part refers 
to adsorption on the reconstructed surface showing (d) the (2x1) missing 
row substrate, (e) (2x2) adsorption (© = 0.25), and (c) (2x1) adsorption 
(0 = 0.5). 

Fig. 3 

Color coded contour plots of the charge density difference, Apva;(r), of Cu(lOO) 
covered with alkah adsorbates (a) Li, (b) Na, and (c) K. The plots refer to the 

(2x1) missing row reconstructed substrate and adsorbate coverages = 0.5 
and show ApvaK^^) ^r a planar cut perpendicular to the surface along the 
substrate troughs. For a definition of Ap^aiir) see text. The units used for 
the color coding, see right scale, are 1/bohr^. 

Fig. 4 

Level diagram of the reconstruction energy Es{Q) for the Cu(100)+Ad, 
Ad = Li, Na, K systems at adsorbate coverages = 0, 0.25, 0.5 (0 = 
denotes the pure substrate). For a definition of Eg see text. The error bar of 
±0.02 eV shown with the Li result for = 0.5 gives an estimate of numerical 
errors in the calculations. 
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